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ABSTRACT 

The near-Earth asteroid (NEA) (175706) 1996 FG3 is a particularly interest- 
ing spacecraft target: a binary asteroid with a low-Af heliocentric orbit. The 
orbit of its satellite has provided valuable information about its mass density 
while its albedo and colors suggest it is primitive or part of the C-complex taxo- 
nomic grouping. We extend the physical characterization of this object with new 
observations of its emission at mid-Infrared (IR) wavelengths and with near- 
IR reflection spectroscopy. We derive an area-equivalent system diameter of 
1.90 ± 0.28 km (corresponding to approximate component diameters of 1.83 km 
and 0.51 km, respectively) and a geometric albedo of 0.039 ± 0.012. (175706) 
1996 FG3 was previously classified as a C-type asteroid, though the combined 
0.4-2.5 yum spectrum with thermal correction indicates classification as B-type; 
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both are consistent with the low measured albedo. Dynamical studies show that 
(175706) 1996 FG3 has most probably originated in the inner main asteroid belt. 
Recent work has suggested the inner Main Belt (142) Polana family as the possi- 
ble origin of another low-Af B-type NEA, (101955) 1999 RQse. A similar origin 
for (175706) 1996 FG3 would require delivery by the overlapping Jupiter 7:2 and 
Mars 5:9 mean motion resonances rather than the uq, and we find this to be a 
low probability, but possible, origin. 

Subject headings: minor planets, asteroids: general 



1. Introduction 

Among the ~ 37 known NEAs with satellites, (175706) 19 96 FG.s has a particularly low 
Av value, making it an ide al target for a spacecraft mission ( jPerozzi et al.ll2001t IChristou 



2003: Johnston et al. 



2OI0I). 



In fa ct, it is the baseline target of the ESA mission study 
Marco Polo-R fiBarucci et al.ll2"oil[ ). (175706) 1996 FG3, hereafter 1996 FG3, is a prototype 
for the "asynchronous" NEA binaries. About 1 5% of all NEAs and small main belt asteroids 
(D < 10 km) are estimated to be such a binary ( iPravec et al.ll2006l ). and are characterized by 
a rapidly rotating primary (P < 4 h) , a nearly spheri cal or oblate primary, and a moderately- 
sized secondary on a close orbit (jPravec et al.ll2006l ). 



1996 FG 3 was the first binary NEA for which eclips e events were detected in optical 
wavel engths (Pravec et al.ll2000l : iMottola fc LahuUal 120001 ). From these and later observa- 
tions, IScheirich fc Praved ( l2009l ) found the period of the mutual orbit to be 16.14 ± 0.01 hr, 
and a diameter ratio of D2/D1 = 0.28!!]o'g2- A circular orbit is consistent with the data, 
but the preferred solution has e = O.lO^gJg and semimajor axis a ~ 3.4 primary radii. The 
primary is an oblate ellipsoid with an axis ratio around a/b ~1.2, while the secondary is 
prolate with an axis ra tio around 1.4. The spin period of the primary is 3.6 hr, that of the 
secondary is unknown. IScheirich fc Praved also find a mass density of 1.4lQ'g g/cm^ and an 
ecliptic latitude of the mutual spin axis of -84°. 

The formation of these binaries points to a history of spin rate increase due to the 
thermal YORP effect, followed by re-shaping and mass- loss ending with a satellite in a 



^Partially based on observations obtained at the European Southern Observatory (ESO, program ID 
383.C-0179A). Observations were also obtained at the Infrared Telescope Facility, which is operated by the 
University of Hawaii under Cooperative Agreement no. NCC 5-538 with the National Aeronautics and Space 
Administration, Science Mission Directorate, Planetary Astronomy Program. 
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close orbit ( jWalsh et al.l (|2008[ ): see also IScheerea (120071 ) work on fission of contact binary 
asteroids). Tfiis process demands a "rubble pile" or gravitational aggregate makeup, where 
bulk re-sha ping is the cause fo r the ubiquitous oblate "top-shape" and equatorial ridge (see 



1999 KW4, lOstro et al.ll2006[ ). The process may end with preferential material migration 



from the poles of the primary toward s the equator and also with general regolith depletion 
(IWalshet aljboosl : IPelbo' et al.lboill ). 



Photometric observations detern iined the optical m agnitude H = 17.76 ± 0.03 and 
phase coefficient G = —0.07 ± 0.02 ( jPravec et al.l 120061 ). B ased on therni a l obse rvations 
(3.6 and 4.5 /xm) using the "Warm Spitzer" space telescope, iMueller et al.l ( 1201 ll ) find an 



area-equivalent diameter of this binary system of l.Slo'g ^^^ ^^^ ^ geometric albedo of 
Pv = 0.041q'q2, consistent with a taxonomic classification in the C complex. It should be 
noted that their diameter and albedo values are based on an assumed color temperature (as 
opposed to a r neasu red one), subjecting their results to significant systematic uncertainty. 
Wolters et al.l (J201ll ). observing at VLT, determined an effective diameter of 1.71±0.07 km 
and a geometric albedo oi pv = 0.044 ±0.004; furth ermore they estimate the thermal inertia 
to be r = 120 ± 50 J m~^ s~°^ K~^. iMueller et al.l also studied the thermal history of 1996 



FG3 over its chaotic orbital evolution finding that it was likely (90%) heated above 482 K, 
possibly hot enough to induce thermal alteration of previously primitive surface material. 



Pravec et al.l (120001 ) and lMottola fc LahuUal (I2OOOI ) estimated that 1996 FG3 is a C-type 



body from its broadband color indices and the steep phase dependence. The second phase 
of the SMASS sur vey's visible wavelength spectroscopy (0.44-0.92 micron) classified 1996 
FG3 as a C-tvpe (IBus Ill999l ) , and when the spectroscopy of 1996 FG3 was extended to 
0.4 - 1.6 /im by iBinzel et al.l ( 1200 ll ) 1996 FG3 rernained classified as a C-type. Recently, 
the taxonor nic classificatio n scheme of iDeMeo et al.l ( l2009l ) h as leveraged data reac hing 2.45 



/im, though iDeMeo et al.l did not initially examine 1996 FG3. iDe Leon et al.l (120111 ) recently 
analyzed three spectra of 1996 FG3, finding multiple acceptable taxonomic classifications; 
Ch, C, Xk and B-type. 

In this study, we present new multi-wavelength observations, refiectance spectroscopy 
between 0.4 and 2.5 /im, and thermal-infrared photometry from 8-19 /im to extend the 
physical characterization of 1996 FG3. We utilize dynamical modeling of main belt source 
regions for NEAs and analyze the possibility of the B-type Polana family as a source for 
1996 FG3. 
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2. Thermal-IR observations and thermal modeling 

We measured the thermal emission of 1996 FG3 from observations at mid-IR wavelengths 
(8-18.5 /im) and also obtained its visible reflected light. Combined with a suitable model of 
the surface thermal emission we determined its size and albedo; results are given in Sect. 12.31 



2.1. IRTF-MIRSI observations 



We observe d 1996 FG.s on 1 M a y 2009 UT using t he Mid-Infrared Spectrometer and 



Imager MIRSI flDeutsch et all l2003l : iKassis et all l2008f ) on the NASA Infrared Telescope 
Facility (IRTF) on Mauna Kea, Hawaii. We used photometric filters centered at 8.7, 9.8, 
11.6, and 18.4 /xm with spectral bandwidths just short of ~ 10%. A standard four-point chop- 
nod pattern was used to remove background flux. The night was clear with n o perceptible 
clouds, and the atmospheric humidity decreased from ~ 30% to ~ 6%. The ICohen et al. 



(I1995I ) flu x calibrator a Hya w as observed before and after the target at very similar airmass. 
A second I Cohen et al.l ( 1l995l ) standard, a Boo, was observed afterwards for cross-checks 
on the flux calibration. The IRTF IDL librarjo was used to coadd data as needed and to 
remove known instrument artifacts. Raw signal counts were derived using standard synthetic 
aperture procedures. Flux calibration factors and atmospheric extinction coefficients for each 
filter were derived from the observati ons of a Hya . Anal yzing the calibrator observations 
(including a Boo) we reproduce their ICohen et al.l ( 119951 ) fiuxes to within 1% for N-band 
fiuxes, and to within 4% for Q-band (18.4 /xm). Color corrections to asteroid fiuxes are 
expected to be at the percent level and were not applied. Final MIRSI fiuxes are given in 
Table [2 

In the N-band filters, we typically performed one or two back-to-back repeat observa- 
tions to ascertain repeatability, which was within 10%. No such variation is seen in repeat 
observations of the fiux calibrators, likely due to fiuctuations in the background level. Our 
observations cover a fraction of the 3.6 hr rotation period (with a low-amplitude lightcurve), 
and given the reported repeatability, no lightcurve correction was attempted. 

MIRSI observations were interspersed with \^-band photomet ry obseryations using the 
optical CCD Apogee, which was calibrated against observations of iLandoltl (119731 ) standard 
SA 101-57. Due to the proximity in airmass (<0.05), extinction correction is not critical; 
we assumed an extinction coefficient of 0.12 mag/airmass. The 12 measured V magnitudes 
are constant within the photometric uncertainty of typically 0.15 mag per data point, no 



^http: //irtfweb. ifa.hawaii.edu/~elv/mirsi_steps . txt 
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indication of a \^-band lightcurve was detected. The average V mag measured is 16.94, in 
excellent agreement with the expected V magnit ude of 16.94 ± 0.13 for lightcurve average 
(based upon H = 17.76 ± 0.03, G = -0.07 ± 0.02, IPravec et al]l2006h . 



2.2. VLT-VISIR observations 



Independent thermal-IR observations of 1996 FG3 were obtain ed on 2 May 2009 UT 
using the VLT Imager and Spectrometer for mid Infrared ( VISIR; iLagage et al.l 120041 ) in- 
stalled at the 8.2m VLT Melipal telescope of the European Southern Observatory (ESO), 
Cerro Paranal, Chile. Photometric observations were carried out through narrow-band fil- 
ters centered at 8.59, 11.88 and 18.72 /im. The obser vation design was largely analogous 
to that of our MIRSI observations (see Sect. 12. ip . The ICohen et al.l ( 1 19951 . Il999l ) standard 
star IID93813 was observed at an airmass very similar to that of the science target and used 
for absolute flux calibration^ Data were reduced as described above, except in the case of 
the PAH_2 filter centered at 11.88 /xm, where coadded asteroid data appeared "smeared." 
In order to correct for the flux lost outside the nominal 10-pixel-radius synthetic aperture, 
the nominal 11.88 /xrn flux was increased by a factor of 13%, derived using the growth curve 



method (JHowelll 119891 ). To test the importance of this approximation we fit the data with- 
out this point, and also with this point with a 10% and 20% uncertainty, finding at most a 
3% change in diameter and no change in rj. See Table [2] for final fluxes. 



2.3. Diameter and albedo results 



We used the Near-Earth Asteroid Thermal Model (NEATM, lHarrislll998l lR to fit the 
thermal fluxes reported in Table [1] and [2l NEATM fluxes are calculated by integrating the 
Planck function over the visible and illuminated part of the surface of a spherical model as- 
teroid. The surface is assumed to be in instantaneous thermal equilibrium with the absorbed 
sunlight, where temperatures across the surface can be rescaled to match the spectral energy 
distribution of the observed fluxes. This temperature rescaling is expressed in terms of a 
dimensionless "beaming parameter" rj such that T^ oc rj. Note that our observations do not 
spatially resolve 1996 FG3. Thus, the NEATM-derived diameter D is that of a sphere with 
the combined cross-sectional area of the two components. The (area-equivalent) component 



see http : //www . eso . org/sci/f acilities/paraxial/instruinents/visir/tools/zerop_cohen_ Jy . txt 



See also iDelbo' fc HarrisI (|2002l ) for an introduction and overview. 
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diameters Di and D2 are related to D via Di'^ + D2^ 



D^ 



In order to determine the statistical uncertainty in the fit parameters of t he NEATM D 
Pv, an d rj, we performed a Monte-Carlo analysis analogous to that described by lMueller et al 
( I2OO7I ) in which 300 random flux sets are generated, normally distributed about the mea- 
sured data. The median of the Monte-Carlo results is adopted as the nominal result with 



asymmetric error bars encompassing the central 68.2% of the results (see lMueller et al.ll2011 
for a more detailed discussion of this method). 

This Monte-Carlo analysis was performed for each dataset (see Tables [1] and [2]) and 
for their combination — note that the observation geometry is essentially identical between 
the two nights; the corresponding change in expected flux is of the order of 1.6%. Given 
the apparent discrepancy between the Q-band fluxes measured during the two nights (see 
below), we repeat this analysis for all datasets minus Q-band fluxes. All results along with 
Monte-Carlo uncertainties are given in Table |3l The data are plotted in Figure [1] along with 
the adopted best-fit model curve. 

The true uncertainties in D and pv are dominated by the systematic uncertainty inher- 
ent in the NEATM, which is conventionally estimated to be 15% in D and 30% in pv (JHarris 
20061 ). In the (exceptional) case of 1996 FG3, the uncertainty in H does not contribute sig- 
nificantly to the error budget. Our results from the different datasets are in excellent mutual 
agreement; in particular, the Q-band fluxes are largely inconsequential for our purposes. We 
adopt the average of the first six lines in Table E] as our final result, including the systematic 
D and pv uncertainty; the adopted 77 uncertainty is the scatter between the six datasets: 
D = 1.90± 0.28 km, pv = 0.039 ± 0.012, and r? = 1 61 ± 0.08. Using the known diameter 
ratio of D2/D1 = 0.281q;q2 (IScheirich fc Praved 120091 ). the corresponding (area-equivalent) 
component diameters are Di = 1.83 ± 0.28 km and D2 = 0.51 ± 0.08 km. 



Our results are in excellent agreement with those reported by lMueller et al.l ( 120111 ) based 



on Warm-Spitzer observ ations (D = 1.8lo;5 ^^ ^^^ Pv = O-^'^-om) ^^^ with those reported 



bvlWolters et aLl teOllh {D = 1.71 ±0.07 km, p^ = 0.044 ±0.004, and r/ = 1.15). Note that 
Wolters et al.r s observations took place at a significantly lower phase angle (11. 7°) than ours 



(67.4 °) , hence the difference in best-fit r] values is consistent with expectations (JDelbo' et al. 
20071 ). iMueller et al.l ( 1201 ll ) did not constrain i], but assumed a value. The good mutual 



diameter agreement o f these three s tudies therefore provides independent support for the 
assumptions made by iMueller et al.l ( 1201 ll ) . 
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3. Visible-Near-IR reflectance and 1996 FGs's main-belt origin 



As part of the MIT-UH-IRTF Joint Campaign for NEO Spectral Reconnaissance tlie 
reflectance of 1996 FG3 was measured between 0.8 and 2.5 /im. In Figure |2] ttie visible 
reflectance from t he SMASS survey is combined with the near-IR to provide reflectance 
from 0.37-2.5 /im (JBinzel et al.ll2004l ). The spectrum presented was measured on 30 March 
2009. A total integration time of 2880 seconds was obtained during a one-hour interval 
beginning at 11:17 UT. A nearby solar analog star, Landolt 105-76, was observed at similar 
airmass immediately following the asteroid measurements. 

The 30 March 2009 spectrum of 1996 FG3 does not yield a unique class, rather the 
allowable classes of C, Ch, or Xk type was determined by t he online visible-near-IR taxonomy 



classi fier based on the DeMeo visible near-IR taxonomy fJDeMeo et al.l 120091 : iDe Leon et al. 



2OIII ). However, a B-type classification in the DeMeo taxonomy demands a negative slope 
to 2.45 /xm, and in the spectrum a thermal tail is evident starting ~2.0 /xm (see Fig. [2]). 
We fit the NEATM to the thermal tail of the SpeX data using the method described in 
Delbo' et al.l (120111 ). A large range in pv^tj combinations fit the thermal tail well, though we 
use py=0.04, rj ~ 1.2. The phase angle during the SpeX observations was ~ 8°, i.e., much 
lower than during our mid-IR photometric observations, hence such a low 77 is consistent 
with our expectations. 

Th e thermally-correcte d spectrum was classifi ed as a B-type by the online classification 
system JPeMeo et al.ll2009[ ). IPe Leon et all J201lh find that the 27 April 2009 spectrum also 
taken by the MIT-UH-IRTF campaign of lower SNR is classified as a B-type by the same 
online tool. However, their spectrum, taken at Telescopio Nazionale Galileo (TNG) on 9 
January 2011, is labeled a Ch or Xk type. We move forward with a B-type classification, 
though naming a taxonomic type is less indicative of composition than an analysis of spectral 
features. Therefore we provide a quantitative comparison and a qualitative interpretation 
of inferred spectral features. 

The ~ 50 km Main Belt asteroid (142) Polana is the largest member of a family p roposed 
to have been the main-belt origin of B-type NEO 1999 RQse (jCampins et al.ll2010l ). due in 
large part to its B-type spec trum. To quant i fy the simi larity between these sp ectra we follow 



the x^ formulation used by iDe Leon et al.l ( 20101) and 



Carapins et al.l (I2OIOI ) between 0.475 



Campins et al.l (120101 ) found a x^ = 4.91 when 



and 1.925 /xm with spacing of 0.0025 /xm. 
comparing 1999 RQse with (142) Polana. Using a similar, but not exact, wavelength range 
we find a value x^ = 4.83 for these two objects. 1996 FG3 compared to 1999 RQse and (142) 
Polana yield values of y^ = 1 .09 an d 1.86, showing a strong similarity to both. The closest 



fit found by ICampins et al.l (120101 ) for 1999 RQse and 27 main belt B-type asteroids was 



X^ 



1.33, so its likeness to another B-type such as 1996 FGs is naturally greater than that 



to the other main beh bodies in that study. Similarly, the match between 1996 FG3 and 
(142) Polana is also closer than for 1999 RQse and (142) Polana. However, when comparing 
spectra via a x^ routine, features representing mineralogical differences can be marginalized 
in favor of bulk shape or slope. 

There are subtle but clearly revealed differences between the three asteroids (Figure [2]) 
that may allow differing mineralogical interpretations. The features near 1.2 and 2.0 /im for 
1996 FG3 are consistent with the presence of some olivine and some pyrox ene, where the 



shallo w depth of the features could be due to opaques. Note, however, that iDe Leon et al. 



( 1201 ll ) rejects these two features due to the shape and location of the fir st band. This 



i s an important disagreement, as our interpretation imply silicates, whereas iDe Leon et al. 



( I2OIII ) does not. The deeper 1.2 /xm feature for Polana, with no feature at 2.0 /xm is 
consistent with the presence of olivine (and perhaps opaques). Having spectral hints of 
olivine and pyroxen e is possibly consistent with the spectral prope r ties o f ureilite meteorites 
(ICloutis et al.ll2010l ). whereas the interpretation of iDe Leon et al.l (1201 ll ) is more suggestive 
of carbonaceous chondrites (they suggest CM2). However, 1999 RQ36 is itself distinct from 
the others for its lack of 1 and 2 /xm spectral features, which is more typical of carbonaceous 
chondrites. 



De Leon et al.l (I2OIII ) also made comparisons with meteorite spectra from the RELAB 



databse. Their TNG spectra was best fit by a combination of CM2 carbonaceous chondrites, 
and L4 and H5 ordinary chondrites. Their preferred match for the high-SNR MIT data was 
with a weakly shocked H4 ordinary chondrite. 



4. Dynamical history and 1996 FGa's main-belt origin 

The chaotic nature of NEA orbital evolution makes it impossible to know the exact his- 
tory for any given body. However, detailed orbital evolution modeling bv lBottke et al.l (120021) 



allows probabilistic estimates for Main Belt source regions to be made. The iBottke et al. 



work finds that the current orbit of 1996 FG3 supports a 92% likelihood that it became an 
NEA after escaping the inner Main Belt via the uq secular resonancqj. 

The z/g marks the inner limit of the main asteroid belt and thus asteroids only reach 
this resonance by migrating inward from orbits with a larger semi- major axis (a > 2.15 AU). 
Similarly, this region is bounded at large semi-major axis by the 3:1 mean motion resonance 



"^ This probability was determined using the 3-source NEO model, where 3 Main Belt sources ar e mod eled 
as the parent populations for NEOs. This is the technique used for 1999 RQ36 bv lCampins et al.l pOlOl ). 
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(MMR) with Jupiter located at 2.5 AU. Thus, escape by the uq resonance requires an initial 
semi-major axis range of 2.15 < a < 2.5 AU, and inward migration. Escape by the z/g 
largely retains the inclination of the orbit so that the current 1.99° inclination is likely near 
its original inclination in the Main Belt. 

A similar analysis has been done for the NEA 1999 RQae where its orbital, physical and 
rotation properties were used to infer a parent family in the inner main belt. Its match with 
the (142) Polana family is due in large part to its classification as a B-type asteroid, and 
similarities in their visible/near-IR spectra. Dynamically, 1999 RQse was found to have a 
95% probability of coming from the inner belt via the uq resonance, and has a low inclination 
of i ~ 6.04°. 1999 RQse is in retrograde rotation, arid also has a "top-shape" which has been 
found frequently among binary NEAs ( iNolan et al. 1 120071 : [Harris et al.ll2009l ). 



Despite uncertainties in the Polana family's age, small members (H > 18.5) have already 
reached the uq resonance by the Yarkovsky effect ( ICampins et al.l 120101 ) . The retrograde 
rotation of 1999 RQse also indicates that its Yarkovsky drift is inward, allowing it to drift 
from the (142) Polana family towards the z/g resonance. 

1996 FG3 however, has H=17.76 and thus should not have had time to reach the z/g 
resonance by the Yar kovsky effect according to this estimate of the Polana family extent 
( ICampins et al.ll2010l ) (Note also that no definitive age is known for Polana due to its mem- 
ber's overlap with the Nysa family). However, it is possible for objects the size of 1996 
FG3 to reach the overlapping Jupiter 7:2 and the Mars 5:9 resonances located at a=2 .2559- 
2.2569 AU and a=2.2542-2.2550 AU jMorbideUi fc NesvornvlE999l : lBottke et alll2007h. This 
region of overlapping resonances is discussed in the Supplementary Material of iBottke et al. 
( l2007l ). but its efficiency as a s ource is uncertain for objects of 2 km such as 1996 FG3. 
Figure S6 of iBottke et al.l ( 120071 ) shows roughly equal numbers of 5 km bodies jumping the 
resonance or getting trapped and excited. 

Smaller objects will have higher Yarkovsky drift rates, likely lowering the fraction that 
get excited onto planet-crossing orbits. We tested the efficiency of resonance-crossing for 
~ 2 km bodies drifting at their maximum Yarkovsky drift rates, which provides a lower 
limit. We used a modified version of swif t_rmv s that applied a con stant change in semi- 
major axis of -10"'^ AU Myr~^ as estimated from lBottke et al.l (120061 ) as an upper limit for 
its Yarkovsky drift rate at this size. We tested a set of 155 massless test particles initially 
located at a Ri2.26 AU. The particles had initial e and i values similar to those of Polana, 
and 21 of the 155 particles had their eccentricities excited to Mars-crossing levels while 
crossing the resonance, and evolving onto Mars-crossing orbits. From here their evolution 
onto NEO orbits will be relatively rapid and leave them virtually indistinguishable from 
asteroids entering NEO-space via the uq resonance. 
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The 21 particles excited out of the inner Main belt all had initial eccentricities above 
0.15, pointing to the resonance's increasing effectiveness for higher e. The Polana family 
itself covers a range from e=0.13-0.17, even though (142) Polana itself has an e=0.136. 
Even with the conservative assumption of the maximum Yarkovsky drift rate for a D=2 km 
asteroid, approximately 10 % of the of family members this size would exit the inner Main 
Belt by the overlapping J7:2 and M5:9 resonances. Therefore, 1996 FG3 could still have 



an ori gin with the Polana family as it fits the other major criteria used by ICampins et al. 



(J2OIOI ) to locate (142) Polana as a possible parent family of 1999 RQse: a B-type asteroid, a 
likely inner-Main Belt origin, low inclination, retrograde rotation. Making this study more 
general, to estimate the number of Polana members of this size in the NEO population, is 
challenging due to the uncertainties surrounding the exact size and age of the Polana family. 



5. Results and Conclusions 

The findings of these observations and dynamical studies are; 

1. 1996 FG3 has a diameter of 1.90±0.28 km and a geometric albedo oipy = 0.039±0.012, 

2. Available data do not uniquely place 1996 FG3 into a single taxonomic class; the best 
data (highest SNR) with a thermal flux correction applied, are most consistent with 
the B-class, 

3. 1996 FG3 is likely ( 92%) to have come from the inner Main Belt via the i^q resonance, 
though it could have escaped by the overlapping J7:2/M5:9 resonances located at 
2.255 AU. 

Together, these results provide a much better picture of the origin of 1996 FG3 and 
raise some tantalizing questions about its relation to the NEA 1999 RQ36, another ideal 
space mission target and primary target of NASA's OSIRIS-REx sample return mission 



flLauretta et al.ll2010[ ) 



Delbo' et al.l ( 120111 ) studied the thermal properties of binary NEAs which have been 
observed in the thermal-IR, finding that statistically elevated rj values, used as a proxy for 
thermal inertia, pointed to cooler surfaces than non-binary NEAs. This was interpreted in 
terms of regolith loss during the binary-formation process. Here, 1996 FG3 was fit with 
an 7] of 1.61 ± 0.08. The large phase angle ~ 68° m akes this less const rai ning, but the 



values are elevated compared to the sample presented in iDelbo' et al.l ( 120111 ) . IWolters et al. 



( 120111 ) obtained r] ~ 1.15 at a phase angle of 11.7°. Taken together, the available thermal-IR 
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data should enable a robust determination of the thermal inertia; that will be discussed, 
however, in a forthcoming paper (Mueller et al., 2012, in preparation). In comparison, 
the thermal inertia of 1999 RQse has been estimated to be around 600 J m~^ s~° ^ K~^ 



(lEmerv et al.ll2010r). B ased on the correlation b etween therrnal ine rtia and object size found 
by lDelbo' et al.l ( 120071 ). and using D ~ 0.58 km ( Nolan et al.ll2007l ). a thermal inertia around 
400 J m~^ s"'' ^ K~^ would be expected, somewhat lower than measured. While 1999 RQae 
is not known to be a binary, its radar shape, oblate with a possible equatorial ridge, is 
suggestive of YORP spinup and reshaping. Like in the case of YORP binaries, this could 



lead to an elevated thermal inertia such as that given by lEmery et al.l (|2010[ ) 



1996 FG3 and 1999 RQse share many physical properties and possibly a similar dynam- 
ical history. As shown in Section [31 there is some uncertainty about the precise taxonomy 
for 1996 FG3, though its spectra is qualitatively similar, with a negative slope from 0.5-2.0 
/im, to that of 1999 RQse- The features for 1996 FG3 at 1.2 and 2.0 /im are an important 
distinction between 1996 FG3 and 1999 RQse- If these features are indicative of silicates 
then 1996 FG3 may have more similarities with ureilites than with carbonaceous chondrites. 
In contrast, the featureless spectra of 1999 RQ36 is suggestive of properties similar to car- 
bonaceous chondrites, with C M typically cite d as cl o se matches for "B- t ypes" depending o n 
the level of thermal processing iDe Leon et al.l ( 120111 ): IClark et al.l (l2010[ ): lHiroi et al.l ( 119961 ). 



It is important to note that I Jenniskens et al.l (|2009[ ) obtained a low signal-to- noise visible 
spectrum of asteroid 2008 TC3 prior to its Earth impact and subsequent recovery as the 
Almahata Sitta meteorite. While their i nitial in t erpre tation of this spectrum suggested F- 
or B-type classification in the system of iTholenI (119891 ). that classification is not confirmed 
with the additio n of n ear-in frared spect r um o f a Almahata Sitta meteorite presented by 
Jenniskens et al.l (12010! ) and iHiroi et al.l (120101 ). While a composite average spectrum of 
different stones of Almahata Sitta reproduces the 2008 TC3 visible spectrum, none of the 
measured lab spe ctra matches t he B- class (or any asteroid class defined over the full spectral 
range defined bv iDeMeo et al.l (120091 )). 



We also note that the F-class itself diverges (fails to exist) for objects who se spectra are 
extended from the visible wavelengths to the near-infrared. As discussed by iDeMeo et al. 
( 20091) . objects that are spectrally similar over visible wavelengths and satisfy the iTholen 
( 19891 ) definition of the "F-class" diverge from one another as their spectra are extended into 



the near-infrared. Thus at present, the visible plus near-infrared spectrum of the Almahata 
Sitta meteorite spectrum fits into no asteroid taxonomy class and closely resembles no other 
currently measured spectrum of any asteroid. 



Jenniskens et al.l (20101) also pursued a parent source family for 2008 TC3 similar to the 
dynamical work of ICampins et al.l (|2010[ ) for 1999 RQ36. The most probable source region 
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was found to be the inner main asteroid belt via t he i/p, resonance ('- ^ 80% p robability): the 



same source region as 1996 FG3 and 1999 RQse- iJenniskens et al.l ( 120101 ) speculated the 
(142) Polana family as a source for 2008 TC3. However the jumbled asteroid- meteorite link 
for the spectrum of 2008 TC3 renders any further connection beyond the scope of this work. 

Thus the implication of this study are that two prime space-mission targets, 1999 RQse 
and 1996 FG3 share similar taxonomies, spectra and possibly even rotation history (YORP- 
spinup). Linking them to the same parent family in the inner-main belt, the Polana family, 
is highly speculative as 1996 FG3 would have had to escape by a minor resonance rather than 
vq. However, it cannot be ruled out entirely as we have shown that the J7:2/M5:9 resonances 
can push asteroids these size onto Mars-crossing orbits and on the path to near-Earth orbits. 
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Fig. 1. — Thernial-IR data of 1996 FG3 (see Tables [T] and [2]) and NEATM model curve using 
the adopted parameters from Table [31 
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Fig. 2. — Visible-near-IR spectral reflectance from 0.4-2.5 /im for 1996 FG3, 1996 FG3 with 
tlie tliermal tail removed, (142) Polana and 1999 RQse (B. Clark et al. 2011 submitted). 
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Fig. 3. — The evolution of massless test particles evolving while a constant change in semi- 
major axis of 10"*^ AU/MYR is applied. The plot shows the eccentricity and semimajor axis 
for each of the 155 particles at the start of the integration (open symbol) and at the end 
of the 100 Myr integration (closed symbols). The vertical dashed lines are the location of 
the Jupiter 7:2 MMR and the Mars 5:9 MMR, while the horizontal dashed line represents a 
Mars-crossing orbit at g = 1.83 AU 
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Table 1: MIRSI observations of 1996 FG3 on 2009 May 1 UT. 1996 FG3 was at a heliocentric 
distance of r = 1.057 AU, observer-centric distance of A = 0.1567 AU, and phase angle of 
a = 67.4°; all values were constant during our observations to within the last quoted digit. 



UT 



Wavelength 



Flux 

Jy 



10 



Flux 



05:51 


11.6 


1.23 ±0.03 


2.75 ±0.07 


05:53 


11.6 


1.29 ±0.03 


2.87 ±0.07 


05:54 


11.6 


1.35 ±0.03 


3.01 ±0.07 


05:59 


8.7 


0.86 ±0.02 


3.40 ±0.08 


06:01 


8.7 


0.77 ±0.02 


3.06 ±0.08 


06:10 


18.4 


1.80 ±0.16 


1.59 ±0.14 


06:17 


9.8 


0.96 ±0.04 


3.00 ±0.13 


06:20 


9.8 


1.05 ±0.05 


3.27 ±0.15 


06:23 


9.8 


0.96 ±0.04 


3.01 ±0.14 


06:25 


11.6 


1.23 ±0.03 


2.75 ±0.06 


06:27 


11.6 


1.21 ±0.03 


2.69 ±0.07 


06:30 


8.7 


0.75 ±0.02 


2.98 ±0.08 


06:33 


8.7 


0.69 ±0.02 


2.74 ±0.09 


06:36 


8.7 


0.76 ±0.02 


3.02 ±0.08 


06:43 


18.4 


1.80 ±0.16 


1.59 ±0.14 


06:48 


11.6 


1.10 ±0.03 


2.44 ±0.07 


06:51 


9.8 


0.90 ±0.05 


2.80 ±0.15 



UT 


Wavelength 


Flux 




Flux 




//m 


Jy 


10 


-i4w/mV^m 


00:56 


8.59 


0.68±0.01 




2.95±0.05 


01:14 


8.59 


0.73±0.01 




2.77±0.05 


00:28 


11.88 


1.19±0.02 




2.53±0.05 


01:05 


18.72 


1.31±0.07 




1.12±0.06 



Table 2: VISIR observations of 1996 FG3 on 2 May 2009 UT. The target was at an observer- 
centric distance of 0.156 AU, a heliocentric distance of 1.053 AU, and a phase angle of 69.1°. 
Note: the 11.88 /im flux quoted here includes the 13% smearing correction discussed in the 
text. 
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Table 3: NEATM fit results for 1996 FG3. Eacli dataset is analyzed with and without Q-band 
data. Adopted final results are given in the last line. The uncertainties in the first six lines 
are purely statistical, whereas the last line includes systematics. 
Data set D{km.) pv 1] 

MIRSI+Q L955E? 0.036+°-°°^ l.GT^Hs 

MIRSI-Q 1.90+^6 0.039+^03 ^-^Kao 

VISIR+Q l.SltH? 0.043^^^^ l-5i;S:J^ 
0-11 n n/1 1+0.003 1 cc+o.u 



VISIR-Q 1.84+;;:^^ OMlt^.Zl l-55_o.o6 

Both+Q 1.941H^ 0.037t°::? 1-671^^ 

Both-Q 1.96;°:°^ 0-036;°:°^^ 1-69;^^ 

Adopted 1.90 ±0.28 0.039 ±0.012 1.61 ±0.08 



